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We analyze the eigenvalues and eigenvectors of the staggered Dirac operator in quenched lattice QCD in the 
vicinity of the deconfinement phase transition using the Liischer-Weisz gauge action. The spectral and localization 
properties of the low-lying eigenmodes show characteristic differences between the Z3 sectors above the critical 
temperature T c . These findings can be interpreted in terms of calorons. 



1. INTRODUCTION 



2. PROPERTIES OF CALORONS 



The chiral phase transition is a main feature of 
QCD. In the instanton picture [|l| the chirally bro- 
ken phase is represented by an ensemble of weakly 
interacting instantons and anti-instantons. Due 
to the fermionic quasi zero modes of the (anti-) 
instantons, a finite density of eigenvalues, pDirac, 
is generated near zero in a (large) volume V, and 
according to the Banks-Casher formula [0 



(qq) = -7rp D irac(0)/V 



(1) 



chiral symmetry is broken. In the chirally sym- 
metric phase instantons and anti-instantons form 
strongly bound molecules, the eigenvalues of the 
Dirac operator move towards the bulk of the spec- 
trum and the eigenvalue density vanishes near the 
origin. Of course one would like to verify the es- 
sential parts of this picture in QCD. Therefore 
we study the localization properties of the eigen- 
modes of the lattice Dirac operator in the neigh- 
borhood of the chiral phase transition. In a pre- 
vious paper || we have investigated this question 
with the Wilson gauge action on smaller lattices. 
Here we present results obtained on larger lattices 
with the Liischer-Weisz gauge action. 
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The temporal extension L t of the lattice is re- 
lated to the temperature T by aT = 1/L t , where 
a is the lattice spacing. Due to the periodic 
boundary condition in the time direction, an in- 
stanton becomes a chain of instantons which one 
calls caloron |l|,Q. The properties of a caloron 
are similar to that of an instanton: a caloron has 
topological charge one and a localized left-handed 
fermionic zero mode. However, the localization 
pattern of the zero mode is different. Whereas in 
the case of an instanton the zero mode is localized 
in space and time, the zero mode of a caloron is 
localized in space but delocalized in time. 

Let us now compare the localization of the zero 
mode in the different sectors of the Polyakov loop 
P. We know that in the deconfined phase P clus- 
ters around the values e l6>p , with 6 P = (real 
sector) and Op = ±2tt/3 (complex sector) [||. Be- 
cause the fermionic action does not share the Z3 
symmetry of the gauge action, the localization of 
the eigenmodes of the Dirac operator can depend 
on the Polyakov loop sectors. Indeed one finds 
that the zero mode ipo of the caloron is more lo- 
calized in the real sector than in the complex sec- 
tor, as can be seen from |||| 

|Vo| 2 (x exp[-2(^ - \0 P \) rT]/r 2 . (2) 

Next we will search for caloron-like configurations 
concentrating on the fermionic properties of the 
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Figure 1. Expectation value of 75 vs. eigenvalue A (left) and inverse participation ratio vs. 1(75)! (right). 
Results for the lowest 10 eigenmodes of « 300 configurations at (3 = 8.45 on a 20 3 x 6 lattice are shown. 



caloron described above. 



4. CALORONS ON THE LATTICE 



3. LATTICE SETUP 

We work in the quenched approximation with 
the Luscher-Weisz gauge action 

S g [U] = Pi J2pl s Re Tr [1 - Upi] + 
Ert 3^ 1* [1 - Urt] + 

PsE P9 l^Tr[l-U pg ], (3) 

where J2 pi , Y^rt anc ^ ^2pg means a summation 
over all plaquettes, 2x1 rectangles and paral- 
lelograms, respectively. The coefficients 02 and 
/3s are computed from (3\ via tadpole improved 
pertubation theory. Further we use the staggered 
Dirac operator 



1 

D = ^2^ a ^ x ) [^x+aU^x) - Sy iX -fiUl(y)] , 

where a fI (x) = (— 1)^1+— +^-1. The eigenvalues 
of D come in pairs of ±iA with A real, so we can 
restrict ourselves to positive A in the following. 
In the continuum limit this action corresponds to 
four quark flavors. From now on we set a to 1. 

We have calculated the lowest eigenvalues and 
eigenvectors for lattice sizes 12 3 x 6, 16 3 x 6 and 
20 3 x 6 below the deconfinement phase transition 
at (3\ = 8.10 (/^Wilson ~ 5.8 0) and above at 
Pi = 8.45 (/^Wilson ^ 6.0 @) with the Arnoldi 
method ||. We will present data for the largest 
lattice at (3 = 8.45. 



When searching for calorons on the lattice with 
the help of their fermionic properties, the zero 
mode of the Dirac operator plays the central 
role. Unfortunately the staggered Dirac opera- 
tor has no exact zero modes, but, as we can see 
clearly in the left plot of Fig. [lj there are for both 
Polyakov loop sectors low-lying eigenmodes with 
a high value of 1(75)! and eigenmodes in the bulk 
of the spectrum with values of (75) around zero. 
We find that the low- lying eigenvalues are quasi 
4-fold degenerate in the sense that the number 
of low- lying eigenmodes is almost always 4, 8, 
12, ... (In Fig. [y this is illustrated by the encir- 
cled eigenmodes, which belong to the same gauge 
field configuration. Note that the total number of 
low-lying eigenmodes of that configuration is 8.) 
So one would expect that those low-lying eigen- 
modes correspond to exact zero modes with defi- 
nite handedness in the continuum. 

Let us now look at the localization properties 
of these approximate zero modes. We define a 
gauge-invariant measure of the localization of a 
quark eigenmode ip^(x) (A is the Dirac eigen- 
value, a a color index), the inverse participation 
ratio 



V 



E^Pa(^) 2 



l2 ' 



(4) 



where V is the number of lattice sites and p\(x 
is the gauge- invariant probability density p\(x) = 




Figure 2. Isosurfaces of p\ (light grey) and p\ (dark grey) for one of the marked eigenmodes of Fig. [j]. 
The left and middle plots show them separately at fixed z, in the right plot they are overlaid at fixed t. 



J2a=i I^a( x )| 2 - For a completely delocalized 
state (all p\(x) the same) one finds I2 = 1, 
whereas a state localized on a single lattice site 
(only one non-zero p\(x)) would have I2 = V. 
In the right part of Fig. [j] we plot I2 vs. |(7s)|. 
We see that the modes in the real sector are, in 
general, more localized than those in the com- 
plex sector. In both Polyakov loop sectors the 
approximate zero modes have a greater value of 
I2 than the modes in the bulk spectrum, which 
means that they are more localized. These are 
exactly the properties we expect for caloron con- 
figurations. 

Finally we consider the localization pattern. 
In Fig. H isosurfaces of p\ and the local chiral- 
ity p 5 x (x) = T,a=i ^x( x Tl^x( x ) are Pitted at 
fixed t and at fixed z for one of the encircled eigen- 
modes. This (typical) mode is clearly localized in 
space but not in time and the isosurfaces for p\ 
and p\ are very similar. This is again as is ex- 
pected for caloron-like configurations. 



5. CONCLUSIONS 

Studying the localization and chirality prop- 
erties of low-lying quark eigenstates in quenched 
lattice QCD, we could characterize semiclassical 
properties of the gauge field configurations with- 
out any cooling. For temperatures above T c we 
found isolated modes with definite handedness. 
They are localized in space but not in time with 



characteristic differences between the Polyakov 
loop sectors. Thus they show essential properties 
of quark states associated with calorons. 
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